Distillation by repeated measurements: continuous spectrum case 
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Repeated measurements on a part of a bipartite system strongly affect the other part not mea- 
sured, whose dynamics is regulated by an effective contracted evolution operator. When the spec- 
trum of this operator is discrete, the latter system is driven into a pure state irrespective of the 
initial state, provided the spectrum satisfies certain conditions. We here show that even in the 
case of continuous spectrum an effective distillation can occur under rather general conditions. We 
confirm it by applying our formalism to a simple model. 
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Distillation procedures aiming at driving quantum sys- 
tems into pure states are relevant tools in the field of 
quantum information and computation [l[. In fact, they 
can be exploited to control the state of a quantum sys- 
tem and play a crucial role to initialize quantum systems. 
Various purification schemes have been proposed Q and 
among them is a state generation strategy based on the 
extraction of a state through repeated measurements Q . 

Indeed, for a generic bipartite quantum system consist- 
ing of two interacting parts X and A, repeated measure- 
ments on one part (X) can strongly affect the dynamics 
of the other (A). In the case of measurements project- 
ing X in its initial state, the dynamics of A is governed 
by an effective evolution operator V T (introduced below) , 
and it has been shown that, when the spectrum of this 
operator is such that its largest (in magnitude) eigen- 
value is unique, discrete, and nondegenerate, then the 
non-measured system A is driven toward a pure state ir- 
respective of its initial condition Q . This distillation pro- 
cedure has been also utilized to produce entangled states 
of a multipartite system A in particular allowing to 
establish entanglement between two spatially separated 
systems via repeated measurements on an entanglement 
mediator 

The above requirements for distillation cannot be satis- 
fied if the spectrum of V T is continuous. Such a situation 
can be found when one or both of the subsystems have a 
continuous spectrum. It has been shown that, although 
the measured part X has a continuous spectrum, V T may 
have a discrete spectrum, provided the spectrum of the 
non-measured part A is discrete j(|. In this paper, we 
investigate the case where V T itself is characterized by a 
continuous spectrum and show that even in such a case 
an effective distillation can happen under rather general 
conditions. 

We prepare X at time t = in a pure state |$), e.g. by 
projecting X to this state by a measurement, while A is 
in an arbitrary mixed state p A (0). The unitary dynamics 
of the total system X+A, governed by the time-evolution 
operator U(t), is interrupted by the measurements per- 
formed on X at intervals r. The measurements are so 



designed to project X onto |<I>). The action of each mea- 
surement is represented by the projection operator 

6 = |*)(*| <g> 1 A . (1) 

The state of the total system after N measurements is 
thus described by 

p T (N) oc [OU(t)] n [\$){$\ <g> p A (0)][U\r)d] N . (2) 

Following we introduce the projected evolution oper- 
ator between two consecutive measurements, 

V T = ($|#(T)|*>, (3) 

so that, after N measurements on X, system A is de- 
scribed by the density matrix 

p A (N) = V T N p A (0)Vj N /P(N), (4) 
P(N)=Ti- A {V T N p A (0)Vj N }, (5) 

where the normalization factor P{N) represents the sur- 
vival probability that X is always found (up to N mea- 
surements) in the state |$) by every measurement and 
thus gives the probability to obtain the state ((3]) . 

Since the operator V T is not Hcrmitian, V T ^ VJ, 
in general, we need to set up both the right- and left- 
eigenvalue problems, V t \ue) = \e\ue) and (ve\V t = 
^e(ve\- Let us assume that the spectrum of the operator 
V T is continuous and nondegenerate, and its eigenvectors 
form a complete orthonormal set in the following sense: 
/ dE\u E )(v E \ = 1 A and (v E \u E >) = S(E - E'). The 
operator V T is expanded in terms of these eigenvectors as 

Vr = j dEX E \u E ){v E \, (6) 
and A is driven by the operator 

t N = J dEX^\u E )(v E \ (7) 
as the measurement is repeated. 
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In order to discuss the possibility to obtain an effective 
distillation of A, we will compare the survival probability 
P(N) and the purity of A, quantified by 



H(TV) = Tr A {p A (N)}, 



(8) 



as functions of the number of measurements TV. We en- 
deavor to find conditions under which an effective distil- 
lation of A, in the sense clarified later, can be achieved 
and consider the cases where the right-eigenvectors are 
orthogonal to each other (ue\ue'} = 5{E — E 1 ) (which 
is actually the case in the examples studied below). In 
such cases, the survival probability JS]) takes the form 



(9) 



P(TV) = J dE\\ E r(uE\PA(0)\uE), 
and the purity of A is expressed as 

n(TV) = J dEdE' \\e\ 2N \X E >\ 2N 

x |(m £ |^(0)|^)| 2 - 



(10) 



Let E* be the value of E at which \Xe\ has its (unique) 
absolute maximum, |A_e|^ =b = 0, and consider the Tay- 
lor expansion of A(E) = — In \\e\ 2 around T5„, |A_g| 2Ar = 

e Nln\X E \ 2 ^ e -NA(E t )-NA"(E,)(E-E t f/2+-_ Notice ^ 

the higher order terms become less relevant as TV in- 
creases, and |Afi| 2Ar becomes well approximated by a 
Gaussian 
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(11) 



Observe here that the Gaussian in (fTTj) becomes narrower 
like A^r oc l/vTV as TV increases, and a narrow band 
around E* is filtered in the spectrum. Putting x — (E — 
E*)/An, the survival probability is shown to behave for 
large TV as 



P(TV) w /(TV) 



dx 



e x ' /2 (u e ,+xA n \pa( 



f(N) ( .9(0) + -g"(0)A 2 N 



\ u E,+xA N ) 

(12) 



where g(y) = (uE,+y\pA{0)\uE,+ y ) , and the purity given 
by (fTU|l becomes 



PW fdxdx' (x , +x ,2 )/2 



n(TV) 



P 2 (TV) 
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(13) 



Here we have introduced h(y,y') = \(uE,+ y \pA(0) x 
\ u E,+y') | 2 and denoted its second partial derivative with 



respect to y as h yy (y,y'). One sees that, for large TV, 
the purity II (TV) approaches 1 with a rate determined by 
A 2 ^ ex 1/TV. System A is thus asymptotically purified 
toward \ue,) by the repeated measurements on X. 

In this process, however, one should pay attention 
to the behavior of the survival probability: the purity 
II(TV) should reach 1 quick enough before the probability 
P(TV) decays out completely. The comparison between 
the rates of the decay of P(N) to and the approach of 
II(TV) to 1 leads to the following optimization criterion in 
order to obtain an efficient distillation of the pure state 
\ue,)- If A(K*) = 0, i.e. |AbJ 2 = 1, the exponential de- 
cay e ~ N ME,) m p(TV) disappears and the decay of P(N) 
is ruled by A^v oc l/VTV. This decay is slower than the 
approach of the purity n(TV) to 1, i.e. A 2 N ex 1/TV. There- 
fore, if the magnitude of the eigenvalue associated to the 
eigenstate \ue,} is \^e,\ = 1, system A is driven toward 
the pure state \ue,) with a rate faster than the decay of 
the survival probability P(N). 

Remark that, when (ue\ue') = &{E — E') is not sat- 
isfied, it can be shown that P(N) and II(TV) at best de- 
crease and increase respectively as A 2 ^ oc 1/TV for large 
TV, so that distillation seems more difficult to achieve. 

Model. We now apply, as an example, the above 
framework to a specific model. We consider a particle 
of mass m interacting with a single cavity mode of fre- 
quency d. This system has been studied in 0] to analyze 
how the repeated measurements on the particle affect 
the dynamics of the cavity mode. In that case, although 
the measured part (particle) has a continuous spectrum, 
the effective dynamics of the cavity mode is described 
by an operator V T characterized by a discrete spectrum. 
Here, the opposite case is investigated, that is, the cav- 
ity mode is repeatedly projected onto its initial state by 
measurements. In this case, as we will show, the effective 
dynamics of the particle is described by an operator V T 
having a continuous spectrum. 

The Hamiltonian describing the system is 



H 



2m 



1 



+ Huj a T a+- )+gp(aJ+a), (14) 



where p is the momentum operator of the particle, d and 
are the annihilation and creation operators of the cav- 
ity mode, respectively, satisfying the commutation rule 
[a, af] = 1, and the real parameter g is the coupling con- 
stant. The dynamics described by the Hamiltonian (fT4")) 
is exactly solvable, and the exact evolution operator at 
time t in the Schrodinger picture is given by [a, 0] 



U(t) 



e -i?xP 2 T /2mh e -iujT(a' i S+l/2) e P(g T cJ -g*a) ^5) 



where £ T = 1 — (2mg 2 / hoj)[l — (sin lot) /lot] and g T = 
g(l - e WT )/hio. 

The projected evolution operator V T defined by © 
strongly depends on the choice of the state |$) on which 
the cavity mode is repeatedly projected at intervals r. 
Nevertheless, it results to be diagonal in the momentum 
representation for any choice of |<fr). It follows that, for 
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any choice of the cavity state |$), the operator V™ has 
a continuous spectrum of the form = J dp Xp \p)(p\, 
where X p is the eigenvalue whose explicit form depends on 
the choice of the cavity state |$) to be measured. In the 
following, we examine how the above general framework 
for distillation with a continuous spectrum of V T works, 
by looking at two cases with different cavity states |$). 
In the first case, |$) is assumed to be a coherent state |a), 
while in the second, a number state \n) . Both cases are of 
interest from an experimental point of view @ . The first 
choice leads directly to a Gaussian form of |A p | 2Ar , while 
in the second case, |A P | 2JV itself is not Gaussian, but it 
becomes well approximated by a Gaussian, and we will 
see that the above general framework actually works. 

As the initial state of the particle, we take a generic 
Gaussian state p~a(0) = J dpdp 1 \p)p pp > (0)(p'| with 



_ (p + p'-2p )2 _ (p-pQ2 > )xo 
e 8(Ap)g S(Ap)gng 11 



xe -M^<P-rt<>**'-*»>\ (16) 



where B = ^{Ax^Ap) 2 ^/ K 2 — 1, p and xo, and 
(Ap) 2 and (Ax)q are the averages and variances of the 
momentum and the position, and IIo the purity of this 
initial state 0. 

Cavity coherent state |$) = \a) (a = |a|e* 7 ). In 
this case, the projected evolution operator reads 



y _ e ~iu)T /2-2i\u\ 2 sin(wT/2)e~ ilJT/2 

x / dp\p)(p\e~ 



(17) 



where b = ~ug* T + e lLJT g T a* . The survival probability 
defined in ([5]) takes the form 



P(N) = 



v/l + 2| 5r | 2 (Ap) 2 7V 

(p -Rc6/|.g T | 2 ) 2 



x exp 



2(Ap) 2 [l + l/2|.g T | 2 (Ap) 2 7V] 



, (18) 



where Keb = — 4(g\a\/hui) sin 2 (wr/2) COS7. In the limit 
of large number of measurements, N l/2|<7 r | 2 (Ap)o, 
P(N) reduces to 



P(N) 



y/2\g r \2(Ap)*N 
1 



x exp 



2(Ap) 2 



Re b 



Po 



(19) 



This formula corresponds to (|12l) . without the second- 
order term proportional to Ajy- , through the relationships 
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Reb 



f(N) = ^ 2 TTA 2 N e- 4Nlal ' 2sin ^ sin2{uT/2 \ 



(p*\pa(0)\p*) = 



-(p,-p ) 2 /2(Ap) 2 



(20) 



One can see that the selected momentum p* , being pro- 
portional to Re b given below (|18|) , can be chosen by tun- 
ing the parameters properly. In particular for 7 = 0, tt, 
the exponential decay in the survival probability, i.e. in 
f{N), can be suppressed, being the final momentum p» 
negative in the first case 7 = and positive in the sec- 
ond 7 = tt. The other parameters can be tuned in order 
to obtain the desired modulus of the final momentum 
p*. Of course, the farther it is from the initial aver- 
age momentum po, the smaller the exponential factor in 
(p*\pa(0)\p*) will be, since it depends on p* — po. 
Let us also look at the evolution of the purity: 



U(N) 



l + 2\g T \ 2 (Ap) 2 N 
l/n 2 + 2|.g T | 2 (Ap) 2 7V 



1 - 



i/ng-i A a 

2(Ap) 2 



(21) 

If the initial state of the particle is pure, IIo = T it 
remains pure after each measurement, while in general, 
the purity of the particle increases as the measurements 
go on. This formula is again consistent with (1131) . 

In Fig. [U the probability P(N) in (03) and the pu- 
rity II(iV) in (f2~T|) are plotted for 7 = tt. The purity 
approaches 1 when the probability is yet far from 0. 

Cavity number state |$) = \n). For n = 1, V T 
results in 

V T =J dp\p)(p\(l-p 2 \g T \ 2 )e~ 3 ^/ 2 

x e - i ^rP 2 T/2mh-p 2 \g T \ 2 /2 ^ (22) 

In this case, our generic formulas in (|11[) suggest 
1 



A 



(p*\Pa(o)\p*) 



1 



e -p 2 /2(Ap) 2 % (23) 
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FIG. 1: (Color online) P(N) (squares) vs. U(N) (diamonds) 
as functions of the number of measurements N, when the 
coherent state \a) of the cavity mode is repeatedly measured. 
The parameters are lot = n/4, \Jrnjhw g — 5, po/Vmhw = 
0.1, (Ap) /Vmho = 0.2, n = 1/V2, and |a| = 1, 7 = tt. 
This plot is independent of xo and (Ax)o- The final selected 
momentum is p» = 2po- 
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FIG. 2: (Color online) Exact (filled marks) and approximate 
(empty marks) behaviors of P(N) (squares) vs. II(iV) (dia- 
monds) as functions of the number of measurements N, when 
a number state |1) of the cavity mode is repeatedly measured. 
The parameters are the same as in Fig. [T] while p t — 0. 



and the asymptotic formula (|12l) for the survival proba- 
bility P(N) for a large number of measurements yields 



P(N) 



1 



,- P 2 /2(A P y 



y/6\g r \*(Ap)*N 



(24) 



On the other hand, (JT3J) gives the asymptotic behavior 
of the purity II(iV) for a large number of measurements, 



U(N) « 1 



i/ng-i A2 



(25) 



which is the same as that obtained in the coherent state 
case in (|2"Tj) . These quantities, the survival probability 
P(N) in dM} and the purity 1T(A) in fl25J|, are plotted 



in Fig. [2] compared with the exact results computed nu- 
merically on the basis of the expression for V T in (|2"2"|) . It 
shows that the exact results and the asymptotic formu- 
las match already after a small number of measurements 
and that the purity n(7V) approaches 1 before the sur- 
vival probability P{N) decays to 0. 

Conclusions. We have studied the distillation pro- 
cess, in which one part of a bipartite system is purified 
by repeatedly projecting the other part onto a certain 
state, in the^case where the dynamics is regulated by 
an operator V T characterized by a continuous spectrum. 
When the maximum of the continuous spectrum \\e\ is 
unique at and the second derivative A"(E*) ^ ex- 
ists, the spectrum |A£| 2Af becomes Gaussian around E* 
as the number of measurements ./V increases, and the 
state is purified to \E 3t ). This purification is optimized, 
if it is allowed to tune the parameters so that |A_g| = 1, 
by which the purity increases faster than the decrease 
of the survival probability. The distillation by repeated 
measurements has been considered to be possible only if 
V T has a discrete spectrum. The present analysis reveals 
that this procedure can be applied to a wider class of 
systems. 
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